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Abstract 

Phase stabilities of Hf-Si-0 and Zr-Si-O have been studied with first-principles and thermodynamic modeling. From the obtained 
thermodynamic descriptions, phase diagrams pertinent to thin film processing were calculated. We found that the relative stability 
of the metal silicates with respect to their binary oxides plays a critical role in silicide formation. It was observed that both the 
HfC>2/Si and Zr02/Si interfaces are stable in a wide temperature range and silicide may form at low temperatures, partially at the 
HfC>2/Si interface. 
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The thickness of Si02 as a gate oxide material in 
advanced complementary metal oxide semiconductor 
(CMOS) integrated circuits has continuously decreased 
and reached the current processing limits [1]. Alternative 
materials with higher dielectric constants, such as HTO2 
and Zr02, are considered as candidates to replace Si02 
for further improvement of their performance [2]. However, 
during the thin film deposition or the subsequent rapid 
thermal annealing, oxides, silicates, and silicides may form 
at the interface since most high-A: materials are metal 
oxides [3, 4]. Among those interfacial phases, silicides are 
detrimental to capacitor performance due to their metallic 
behavior[5]. In this regard, thermodynamic stability cal- 
culations and experimental results have shown that the 
interface between HTO2 and Si is found to be stable with 
respect to the formation of silicides [4]. On the other hand, 
the Zr02/Si interface was found to be unstable around 
1000K, which is in contradiction to the calculation by 
Hubbard and Schlom [2]. It was also observed that the 
Hf-silicide forms upon decomposition of Hf02 in low oxy- 
gen partial pressures [5, 6, 7, 8] and HfSi04 suppresses 
Hf-silicide formation [9]. 

Although the phase stabilities in the Hf-Si-0 and Zr-Si- 
systems are important, comprehensive thermodynamic 
explanations are not yet available. In this paper, based on 
the recently developed thermodynamic descriptions of the 
Hf-Si-O[10] and Zr-Si-0 systems with first- principles cal- 
culations and thermodynamic CALculation of PHAsc Dia- 



grams (CALPHAD)modeling[ll], various phase diagrams 
pertinent to thin film processing are investigated. 

In the CALPHAD approach, the Gibbs energies of in- 
dividual phases in a system are evaluated from the exist- 
ing experimental data with the so-called sublattice model 
based on the crystal structures. The Gibbs energies of a 
higher-order system can be readily extrapolated from the 
lower-order systems, and any new phases of the higher- 
order system can be introduced. However, it is not always 
possible to have enough experimental data for thermody- 
namic modeling of a system [2] so that theoretical calcu- 
lations, such as first-principles calculation results, can be 
used as supplementary experimental data. The Hf-Si-0 sys- 
tem was recently modeled with first-principles calculations 
and the CALPHAD approach[10]. The formation enthalpy 
for HfSi04 is calculated from first-principles calculations 
since no experimental measurement is reported. The ref- 
erence states of the formation enthalpy for HfSi04 are de- 
rived from the two binary metal oxides as shown in Eqn. 1, 
where E represents the total energy of each phase. The for- 
mation entropy of HfSi04 was evaluated from the temper- 
ature of peritectic reaction, Hf02 + Liquid — > HfSi04, in 
the Hf02-Si02 pseudo-binary system. The thermodynamic 
description of the Zr-Si-0 system was obtained by combin- 
ing the previous modelings[12, 13, 14] and first-principles 
calculation of ZrSi04 in the present work. 



Aif HfSi0 4 
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£(HfSi0 4 ) - \E{m0 2 ) - i£(Si0 2 ) (1) 

The highly efficient Vienna A b initio Simulation Package 
(VASP)[15] was used to perform the density functional the- 
ory (DFT) electronic structure calculations. The projector 
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augmented wave (PAW) method[16] was chosen, and the 
generalized gradient approximation (GGA)[17] was used to 
take into account exchange and correlation contributions. 
An energy cutoff was constantly set as 500 eV for all the 
structures, and the Monkhorst-Pack scheme was used for 
the Brillouin-zonc integrations. For the fc-point sampling, 
authors aimed all the structures to have the fc-point meshes 
as close as (# of atoms in a structure) x k x x k y x k z ~ 
5000 fc-points. Thus, HfSi04 and ZrSiO<4, for example, have 
8x8x8 fc-point meshes. The calculated results of metal 
oxides and silicates are listed in Table 1. 

From the constructed thermodynamic databases of the 
Hf-Si-0 and Zr-Si-0 systems, the isopleths of Hf02-Si and 
ZrOvSi are calculated in order to investigate the stability 
range of silicides at the metal oxides/silicon interface and 
arc given in Figure 1. Calculated results show that HfSi2 
is stable up to 544K based on the formation enthalpy of 
HfSi04 from first-principles calculations. It is generally ac- 
cepted that the uncertainty of the formation enthalpy of in- 
termetallic compounds, which originates from the density 
functional theory itself, is about ±1 kJ/mol-atom[18, 19]. 
Thus the associated decomposition temperature of HfSi2 
at the Hf02/Si interface ranges from 382 to 670K when the 
formation enthalpy of HfSi04 is adjusted within its uncer- 
tainty range from -0.769 to -2.769 kJ/mol-atom. The for- 
mation entropy of HfSi04 with respect to the binary oxides 
was evaluated correspondingly to reproduce its pcritectic 
reaction at 2023K. It should be noted that the phase sta- 
bility range of HfSi2 in the Hf02-Si isopleth is not directly 
correlated with the first-principles calculation of HfSi04, 
but predicted from the Gibbs energies of other phases, in- 
cluding the HfSi04 phase, in the Hf-Si-0 system. Even with 
the uncertainty of formation enthalpy for HfSiC>4, the tem- 
perature range for the Hf02 and Si coexistent phase region 
in the isopleth is fairly wide from 670 to 1700K. 

For ZrSi04, besides the uncertainty of formation en- 
thalpy from first-principles in the present work, the pcri- 
tectic reaction (Zr02 + Liquid — > ZrSiO^) temperature in 
the Zr02-Si02 pseudo-binary is also uncertain from 1910 
to 1949K. Thus, the formation entropy of ZrSi04 varies ac- 
cordingly. The Gibbs energy of ZrSi04 at 1000K evaluated 
from the formation enthalpy derived from first-principles 
and formation entropy evaluated from the peritectic tem- 
perature of 1949K (listed in Table 1) is almost identical 
to the value used by Hubbard and Schlom [2] . With these 
formation enthalpy and entropy values of ZrSi04, ZrSi2 is 
completely suppressed by ZrSi04 and does not show up in 
the Zr02-Si isopleth. To make ZrSi2 appear in the Zr02-Si 
isopleth, the formation enthalpy of ZrSi04 should be more 
negative than the first-principles calculation result within 
the uncertainty of formation enthalpy and pcritectic tem- 
perature for ZrSi04. When formation enthalpy of ZrSi04 
with respect to the binary metal oxides is set to its lowest 
limit from the uncertainty of first-principles calculations, 
3.358 kJ/mol-atom, and entropy of forma- 

0.788 J/mol-atom-K, ZrSi 2 



AH 



ZrSi0 4 
/ 



is stable up to 879K in the ZrG^-Si isopleth. Then forma- 
tion enthalpy of ZrSi0 4 is AHf si ° 4 = -338.568 kJ/mol- 
atom with respect to SER (Standard Element Reference) 
and this agrees well with the experimental measurement, 
-339.033 kJ/mol-atom from Ellison and Navrotsky [20]. 
Consequently, the safe temperature range for Zr02 to be 
stable with Si is between 879 and 1630K, narrower than that 
of Hf02 and Si. However, even with these uncertainties, 
both metal oxides are stable with Si approximately above 
900K as summarized by Hubbard and Schlom [2](1000K). 
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Fig. 1. Calculated isopleths of Hf02-Si and Zr02-Si at 1 atm. Poly- 
morphs of metal oxides for Hf02 and Zr02, i.e. monoclinic, tetrag- 
onal, and cubic, are given in parentheses. 

These isopleth calculation results are in agreements with 
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Table 1 

First-principles calculation results of metal oxides and metal silicates. 

Phases Space Lattice parameters Total energy Aiij a ASf a ' b 



Group a b c a (3 ■y (eV/atom) kj/mol-atom J/mol-atomK 



Hf0 2 P2 x /c 5.135 5.194 5.314 90 99.56 


90 


-10.2101 






Zr0 2 P2 x /c 5.221 5.287 5.398 90 99.63 


90 


-9.5376 






Si0 2 P3 2 21 5.007 5.007 5.496 90 90 


120 


-7.9581 






HfSi0 4 /4i/amd 6.616 6.616 6.004 90 90 


90 


-9.1024 


-1.769 


-0.219 


ZrSiQ 4 J4i /amd 6.698 6.698 6.038 90 90 


90 


-8.7723 


-2.358 


-0.275 



a Formation enthalpies and entropies of metal silicates are expressed with respect to their binary oxides. 
b Formation entropies are evaluated from temperature of peritectic reactions (M0 2 + Liquid — > AfSiCU). 



Gutowski ct al. [4] for the Hf-Si-0 system but not with the 
Zr-Si-0 system. In their calculations, they assumed that 
thermal effects arc of secondary importance for Gibbs en- 
ergy change so that the contribution from entropy was ig- 
nored in the calculations for silicide formation reactions. 
However, our calculations from the individual thermody- 
namic databases showed that such an entropy effect can- 
not be neglected. According to our calculation results, both 
Hf02 and Zr02 are stable with Si and this is in agreement 
with the calculation from Hubbard and Schlom [2]. How- 
ever, ZrSi2 was found at the metal oxide/Si interface in 
their experiment [4] while Hf02 was stable on Si without 
any silicides formation when they deposited at 823K and 
then annealed at 1073K. It can be explained that since their 
fabrication process was rapid thermal chemical-vapor de- 
position (RTCVD), it might not have reached the thermo- 
dynamic equilibrium state. Furthermore, the oxygen par- 
tial pressure of their experiment was not reported. The ef- 
fect of oxygen partial pressure will be discussed later in this 
paper. 

The calculated isoplcths indicate that metal silicates play 
an important role in the silicide formation as suggested by 
Takahashi ct al. [9] From first-principles calculations, the 
formation enthalpy of ZrSi04 is —2.358 ± 1 kJ/mol-atom 
whereas that of HfSi04 is only —1.769 ± 1 kJ/mol-atom 
when the reference states are set to the binary oxides. It is 
intriguing to see that such a small (0.6 kJ/mol-atom) differ- 
ence in the formation of metal silicates greatly affects the 
phase stability at the metal oxides/silicon interface. This 
can be explained by comparing the relationship between 
metal oxides, silicates, and silicides in the isothermal sec- 
tion. 

Isothermal sections of the Zr-Si-0 system at two dif- 
ferent temperatures, 500K and 1000K, are calculated (see 
Figure 2) to investigate the phase relationship regarding 
the decomposition of ZrSi2 at the Zr02/Si interface. The 
two different three-phase regions, ZrSi04+Zr02+ZrSi2 and 
ZrSi04+ZrSi2+Si, in the 500K isothermal section are in- 
tersected by the line connecting Zr02 and Si. Therefore, 
ZrSi2 can be found in the thin film process. However, the 
1000K calculation shows that Zr02 is stable with the Si 
substrate without any silicide formation as there is a tie line 
connecting Zr02 and Si. Isothermal sections of Hf-Si-0 at 
the same temperatures, 500K and 1000K, showed similar 
phase stabilities as Zr-Si-0 [10]. 
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(b) Zr-Si-O at 1000K 

Fig. 2. Calculated isothermal sections of Zr-Si-O at (a) 500K and 
(b) 1000K. Tie lines are drawn inside the two phase regions. 

According to the calculation results of isoplcths and 
isothermal sections of the Hf-Si-0 system, Hf02 and Si 
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should be stable at the temperature range between 670K 
and 1700K. However, it is reported that under oxygen- 
deficient conditions, Hf-silicidc forms at the HfC^/Si inter- 
face even in this temperature range. Wang et al. [6] found 
that oxygen-deficient HfO K <2 consumes the oxygen in the 
SiC>2 thin layer covered on the silicon substrate to form 
fully oxidized metal oxide. Even the layer of silicates will 
be decomposed along with the progress of HfO x <2 deposi- 
tion. The recent work from Miyata et al. [7] also confirmed 
the formation of nanometer-scale HfSi2 dots on the newly 
opened void surface produced by the decomposition of 
Hf02/SiC>2 films at the oxide/void boundary in vacuum. 

To further understand the effect of the oxygen partial 
pressure, the phase diagrams of oxygen partial pressure- 
temperature are calculated and shown in Figure 3 with the 
ratio between the metals (Hf and Zr) and Si set to 1. It 
should be mentioned here that the oxygen partial pressure 
in these calculations are the local oxygen pressure at the 
interface, which is extremely low. Consequently, both sys- 
tems initially form only metallic silicides. As oxygen par- 
tial pressure increases, part of the silicides transform into 
metal oxides. Afterwards, HfC>2 and Si arc stable in the Hf- 
Si-0 system as confirmed by experiments. [5, 6, 7, 8] Then, 
HfC>2 is in equilibrium with HfSiO.4. In the Zr-Si-0 system, 
with further oxidization, ZrC>2 is stable with Si. However, 
the phase region is narrower than that of the Hf-Si-0 sys- 
tem. This is in agreement with Copcl et al. [3] that ZrC>2 is 
vulnerable to high temperature vacuum annealing. There- 
fore, it is possible to have a stable ZrO^/Si interface, but 
this is very challenging in the high vacuum condition. 

In summary, with the thermodynamic descriptions of the 
Hf-Si-0 and Zr-Si-0 systems developed by the CALPHAD 
technique, isoplcths and isothermal sections can be readily 
calculated. It is found that the HfC^/Si interface is thermo- 
dynamically stable between 670 and 1700K as far as oxy- 
gen partial pressure is high enough to keep Hf02 stable. 
Zr02/Si interface is stable between 879 and 1630K, but in 
the oxygen-deficient condition, the processing window for 
a stable ZrO^/Si interface is very narrow. Both metal ox- 
ides are stable with a Si substrate above 900K, even with 
the uncertainties of the formation enthalpies and entropies 
for metal silicates. 
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